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It is analyzed whether the potential energy landscape of a glass-forming system can be effectively 
mapped on a random model which is described in statistical terms. For this purpose we generalize 
the simple trap model of Bouchaud and coworkers by dividing the total system into M weakly 
interacting identical subsystems, each being described in terms of a trap model. The distribution 
of traps in this extended trap model (ETM) is fully determined by the thermodynamics of the 
glass-former. The dynamics is described by two adjustable parameters, one characterizing the 
common energy level of the barriers, the other the strength of the interaction. The comparison is 
performed for the standard binary mixture Lennard- Jones system with 65 particles. The metabasins, 
identified in our previous work, are chosen as traps. Comparing molecular dynamics simulations of 
the Lennard- Jones system with Monte Carlo calculations of the ETM allows one to determine the 
adjustable parameters. Analysis of the first moment of the waiting distribution yields an optimum 
agreement when choosing M ~ 3 subsystems. Comparison with the second moment of the waiting 
time distribution, reflecting dynamic heterogeneities, indicates that the sizes of the subsystems may 
fluctuate. 



PACS numbers: 0.8.15 

I. INTRODUCTION 

Understanding the properties of supercooled liquids 
and in particular the nature of the glass transition is 
still a major scientific challenge. An important concept 
to grasp the important physics is related to the poten- 
tial energy landscape (PEL) HHQ. In this view the 
system is regarded as a point in the high-dimensional 
configuration space. At sufficiently low temperatures the 
system is mainly characterized by the different minima 
{inherent structures) of the PEL and their respective har- 
monic attraction basins 

H H H With the advent 
of faster computers in recent years it became possible 
to elucidate the properties of the minima and the sad- 
dles and thus to relate the thermodynamic 0, 0, IS [H3| 
as well as the dynamic behavior of supercooled liquids 

EH E3, E3, E3, EE EE E3, E3, ISH H3, IS3i to the 

properties of the PEL. Most of the work has been devoted 
to the study of the thermodynamics. From analyzing the 
PEL for different densities it was possible, e.g., to ex- 
tract the equation of state over a very broad range of 
temperatures and pressures 0, 13 . 

The Adam-Gibbs relation constitutes a relation be- 
tween the dynamics and the thermodynamics j2f|. It 
expresses the relaxation time in terms of the configura- 
tional entropy. This relation seems to hold rather well for 
different systems 0, |2(| H3 • The theoretical relevance, 
however, is not clear so far. Therefore it is still an open 
question which properties of the PEL beyond the dis- 
tribution of inherent structures determine the dynamic 
behavior |28j . 



It has been attempted to fill this gap between ther- 
modynamics and dynamics by phenomenological mod- 
els for which the system is hopping between discrete 
states like for the random-energy model or the trap model 
HE HE EIIl • Here we are in particular interested in the 
trap model as discussed by Monthus and Bouchaud |30| . 
This model starts from a distribution G(e) of traps with 
depth e. The average waiting time in a trap of depth e is 
given by r(e, T) = r exp(/3V(e)) where V{e) — b — e and 
b is the common barrier level. After thermal excitation 
the system randomly chooses a neighbor trap; see FigQ] 
for a one-dimensional sketch. The trap energies are ran- 
domly distributed. Of course, for a real system one has 
to consider a distribution of traps in high-dimensional 
space. At low temperatures the slow dynamics is related 
to very long residence times in deep traps. Most work 
has been devoted to exponentially distributed traps but 
different distributions G(e) can be chosen as well. In any 
event, G(e) fully determines the thermodynamic proper- 
ties and should reflect the energy-distribution of inher- 
ent structures. A priori, it is not clear whether the trap 
model is of major relevance for describing the PEL of 
real glass- forming liquids. From a conceptual point of 
view an important property of the trap model is its non- 
topographic nature. The arrangement of traps is purely 
statistical and the dynamics is exclusively related to the 
properties of the individual traps and not to possibly 
subtle correlations among adjacent traps. We mention 
in passing, that in alternative approaches the dynamics 
is rationalized without any reference to thermodynamic 
properties }32| . 
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configurational coordinate 



Figure 1: Thick line: Sketch of the standard trap model. 
Thin line: The traps are identified as the metabasins and 
typically contain several inherent structures. 

A very detailed analysis of the dynamics naturally in- 
volves the properties of the barriers between the inher- 
ent structures as well as the topology of the PEL. For a 
binary-mixture Lennard Jones system (BMLJ) trajecto- 
ries, generated via molecular dynamics simulations, have 
been analyzed in detail. It is immediately visible that 
during very long times the system resides in some fixed 
region of the configurational space. Thus the dynam- 
ics is restricted to simple back and forth jumps between 
a finite number of inherent structures [T^, l2ll l33f | . It is 
useful to regard such a set of inherent structures as an el- 
ementary state, denoted metabasin (MB) [20II34I I35ll36| . 
The energy and configuration of a MB is taken from the 
inherent structure of this MB with the lowest energy. Fi- 
nally, the total trajectory of the system in configuration 
space can be regarded as a sequence of different MB, 
each characterized by an energy and a waiting time. Ac- 
tually, some work has been recently performed to relate 
this configurational space picture to real space dynamics 
[ill Ell HE HI! Oil an d to explain the response to applied 
strain in glasses with the help of the concept of MB [4L| . 

Here we would like to remind the reader of a remark- 
able result, presented in previous work |2l). It turns out 
that for all relevant temperatures analyzed so far, 

N 

Y,((ri(n)-n(0)) 2 )*a 2 n. (1) 
i=i 

The left term denotes the mean square displacement after 
n MB-transitions and a 2 is the temperature-independent 
average distance between two adjacent MB in configura- 
tion space. Two important consequences follow from Eq. 
[TJ (i) On the level of MB the dynamics of the BMLJ sys- 
tem can be described as a random walk in configuration 
space; (ii) The diffusion constant can be expressed as [2l| 

D(T)= a y6(T(T)). (2) 

Thus knowledge of the average waiting time is sufficient 
to predict the macroscopic dynamics. Actually, the sim- 
plicity of these results can serve as the a posteriori justi- 
fication for the introduction of the MB approach because 
on the level of inherent structures no simple physical pic- 
ture can be formulated for the dynamics of the BMLJ 



system |2fj. Recent work for a pure Lennard- Jones sys- 
tem of significantly smaller size seems to indicate that 
there the random-walk like dynamics already holds on 
the level of inherent structures 01 • 

Another remarkable result is the fact that a BMLJ 
system with N = 65 particles has no relevant finite size 
effects for the dynamics in the accessible temperature 
range for computer simulations |43j . In particular, a sys- 
tem with N — 130 particles behaves like an indepen- 
dent superposition of two systems with N — 65 particles. 
Therefore it is sufficient to study the PEL of the BMLJ 
system with A?" = 65 particles. Actually when choosing 
even smaller systems (e.g. N = 40) significant finite-size 
effects occur |g. 

Maybe the simplest disordered model, reproducing the 
observed features, is the trap model as introduced above. 
Since the distribution of traps is fully determined by the 
thermodynamics of the system, the dynamics is basically 
governed by a single parameter, namely the barrier height 
b. In any event, if the PEL can be mapped on a trap 
model, the MB have to be identified as traps; see Fig^ 
Recent simulation work for the BMLJ system directly 
shows in different ways that on a quantitative level the 
predictions of the trap model are neither compatible with 
the distribution of waiting times |19| nor with the acti- 
vation energy V(e) = b — e ppf . 

The scope of this work is twofold. First, we will use 
general arguments to define an extended version of the 
trap model. It fulfills some basic requirements which any 
model for glass-forming systems should fulfill. Second, 
by detailed comparison with the properties of the BMLJ 
system we elucidate to which level this extended trap 
model will be able to explain the BMLJ dynamics. From 
the remaining differences important information about 
the nature of the PEL can be derived. Among other 
things the concept of dynamic heterogeneities is related 
to properties of the PEL in a detailed way. Finally, some 
possible applications of the present results will be indi- 
cated. 



II. THE BMLJ SYSTEM 

A. Definition 

The BMLJ system is one of the standard glass-forming 
systems used in computer simulations. It contains 80% 
large and 20% small particles. We use the potential pa- 
rameters as outlined in Ref. |20l |44| . We have performed 
molecular dynamic simulations for temperatures between 
T = 0.4 and T = 1 with system size N = 65 and us- 
ing periodic boundary conditions. Here we present data 
for T > 0.45. All units are given in LJ units. The 
mode-coupling temperature for this model is T c « 0.45 
[2ll l4l| . As shown in previous work this system size is 
large enough to reproduce the macroscopic behavior of 
the BMLJ at least for T > T c . 
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B. Some previous results 

The most important quantity for the thermodynamics 
is the effective density of MB, i.e. G(e). In previous work 
it turned out that to a very good approximation 

G(e) cxexp(-( e -e ) 2 /2^o) (3) 

with en = —279.2 and Cq = 9-05 in Lennard- Jones units 
|oll(:ill2(|. From entropy arguments one can estimate that 
G(e) has a low-energy cutoff around e = —306. From 
G(e) one can calculate the Boltzmann probability p(e, T) 
that at some randomly chosen time the system visits a 
MB of energy e at temperature T via 

p(e,T)cxG(e)exp(-e/T). (4) 

Note that G(e) is slightly different from the real density 
of MB due a small dependence of the local curvature 
of the different MB on their energy e [5J, |6j. We just 
mention in passing that basically the same results for 
the thermodynamics are obtained if inherent structures 
rather than MB are analyzed. 

In previous simulations we have analyzed in detail the 
waiting time of the system in MB of energy e. Averag- 
ing over different MB with energy e we have found the 
relation 

(r(e,T))=ro(e)exp(/3^ // (e)), (5) 

i.e. the average waiting time in a MB of energy e shows 
simple Arrhenius behavior. The function V e //(e) can be 
interpreted as the average saddle height when leaving a 
MB of energy e. This function is reproduced in Fig0 
Note that it was possible to identify the saddles in the 
PEL which give rise to the specific values of V e ff (e). For 
e < —295 the function V e //(e) is linear in e with slope 
-0.55. In the high-energy limit V e ff{e) approaches a con- 
stant close to 1. This reflects the fact that the dynamics 
in high-energy states resembles a non-activated fluid-like 
dynamics. Using this choice of V e //(e) it is possible to 
formally describe this fluid-like dynamics as transitions 
between MB |20j . A similar limiting value is reported 
in 46]. In the same energy range the effective prefactor 
To (e) increases by one order of magnitude [2(j ■ 

As shown in |22l | the dynamics is governed by activated 
processes for T < 0.6. Correspondingly, the average wait- 
ing time is dominated by the escape from low-energy MB 
(e < —295). For temperatures below T c it becomes diffi- 
cult to obtain equilibrium data. Therefore we restrict our 
detailed analysis to temperatures between 0.45 and 0.6 to 
compare the dynamics of the BMLJ system with the pre- 
dictions of the modelling. In Fig [21 we have also included 
the distributions p(e,T) in this temperature range, ob- 
tained from EqQ] On the low-energy side the contri- 
bution stems from few events with relatively large wait- 
ing times. Thus energy-dependent observables for, say, 
e < —300 possess rather large statistical errors. In the 
detailed analysis of the subsequent sections we consider 
MB with e 6 [-300, -295]. 




Figure 2: The function V e / / (e) as obtained for the BMLJ sys- 
tem together with the distribution p(e, T) for three different 
temperatures. 



III. EXTENDED TRAP MODEL: DEFINITION 
AND INTERPRETATION 

A. General aspects of statistical models 

How do properties of the PEL translate into dynam- 
ics? According to the landscape paradigm the properties 
of glass-forming systems are reflected by the properties 
of the minima. As discussed above, we would use the 
MB rather than the minima as the basic objects. In any 
event, one has an exponentially large but finite number of 
states. In the most detailed description the dynamics out 
of some state i is governed by two temperature dependent 
functions. First, the probability function r(...,i — > j) de- 
scribes the probability that the next state visited is state 
j. The dots indicate that in principle this probability 
may also depend on the history before entering state i. 
Second, the waiting time in state i is characterized by the 
distribution function of the waiting times /(r|...,i — > j) 
when going from state i to state j. 

Of particular interest are the low-energy states because 
they determine the dynamics at low temperatures. As 
shown in [23 the escape out of these states is activated. 
This implies that the dynamics can be very well described 
as Markovian, i.e. the escape out of state i does not de- 
pend on the previous states. Furthermore the waiting 
time distribution is close to exponential, i.e. character- 
ized by a single average waiting time. This is explicitly 
shown in Fig^for two randomly chosen MB, both having 
energies close to e = —302. To obtain one set of data, we 
have performed repeated simulations from the same ini- 
tial configuration, using different initial velocities. Per- 
forming our standard MB analysis we have identified for 
every run the time when the system is leaving the MB. 

As a consequence, the dynamics is fully characterized 
by the probability function r(i — > j) and the average 
waiting time of state i. Alternatively, the information 
can be expressed by the transition rates T(i — > j). On 
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Figure 3: The probability distribution of waiting times for 
two randomly chosen MB. The time-axis has been normalized 
by the respective average value T averg . A perfect exponential 
waiting time distribution with slope -1 has been included. For 
better visibility the curves are shifted against each other. 60 
escape events have been analyzed for each MB. 



this level, the total dynamics can be formulated as a set of 
rate equations. The rates can in principle be determined 
from simulations. Actually, for very small systems where 
most relevant states and corresponding saddles can be 
explicitl y d etermined, rate equations have been already 
studied |4jj. We mention in passing that in this frame- 
work equilibrium as well as non-equilibrium phenomena 
are fully characterized. 

This approach has to fail for an exponentially large 
number of relevant minima. Rather one has to resort to 
some statistical approach. This involves two steps. First, 
one has to characterize each state i by some relevant in- 
ternal parameters Xi, 1,2^,2, ■■• • As a natural choice en- 
ergy may be taken as one of these parameters. Then a di- 
rect connection to the thermodynamics becomes possible. 
Second, one has to define the transition rates r<r(i — > j). 
In the trap model, e is the only internal parameter and 
Tt(« — ► j) does not depend on j. For a good statisti- 
cal model it would be possible to reproduce any time- 
dependence of the energy for any temperature schedule, 
i.e. also in non-equilibrium situations. Reproducing ob- 
servables like the alpha-relaxation time or the diffusion 
constant requires additional information about the real- 
space nature of the different states. Fortunately, for the 
BMLJ systcm(and also for BKS-SiC>2 (unpublished re- 
sults)), Eq. n expresses a very simple relation between 
the configuration space dynamics and the temperature 
dependence of the diffusion constant. 

Searching for a very simple statistical model one may 
naturally end up with a trap model because it only in- 
volves the energy as the only internal parameter. Be- 
yond its simplicity it has the advantage that of automat- 
ically reproducing the random-walk like nature of dy- 
namics (Eq^). However, as mentioned in the Introduc- 




Figure 4: (a) The division of the total system in M inde- 
pendent subsystems (here: M=3). Of course, the borders 
are expected to fluctuate with time, (b) Every subsystem 
is described by a trap-model, which is characterized by the 
function V(e). 



tion major discrepancies between the BMLJ data and 
the predictions of the trap model have been reported. In 
the following we will formulate an extension of the trap 
model by including more than one internal parameter. 
Then it will be possible to reproduce several non-trivial 
aspects of the BMLJ dynamics. 



B. Relaxation in a collection of elementary trap 
systems 

It is well accepted that for supercooled liquids there 
exists a finite dynamic correlation length, sometimes de- 
noted length scale of dynamic heterogeneities EME^,H3- 
It can even be measured experimentally |5lL |52( • At a 
given temperature, on average P particles (P may be 
somewhat temperature dependent) relax in a collective 
fashion. Stated differently, relaxation processes are re- 
stricted to small regions of the system, containing on av- 
erage P particles. Depending on the exact definition of 
the length scale of dynamic heterogeneities one expects 
values between 10 and 100 particles. Thus one may be 
tempted to divide the total system into M = N/P sub- 
systems which to first approximation behave indepen- 
dently Of course, the division into M subsystems 
should not be taken too literally. Rather these border- 
lines may fluctuate with time (see also below). Guided 
by the random- walk nature of the dynamics in configura- 
tion space we approximate every individual subsystem by 
a trap model; see FigQ]}. The resulting model is termed 
extended trap model (ETM) . Then the total energy e may 

be written as e = X)m=i e ™- A simple trap model would 
correspond to M = 1. However, in this case one would 
have expected that V e //(e) has a slope of -1 for low e 
in contrast to the real behavior; see FigI3 Thus it is 
already clear that a mapping on a simple trap model is 
not possible and that, if at all, M > 1 has to be chosen 
for a BMLJ system with N=65 particles pH l2fj||. 

In order to reproduce the thermodynamics of the 
BMLJ system the distribution of traps G(e) in each sub- 
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system has to be chosen as 

/ (e m -(e /M)H , fi v 
GM(e m )ocexp^ 2(g 2 /M) j • (6) 

The dynamics is expected to strongly depend on the 
length scale of dynamic heterogeneities, i.e. on the value 
of M. In what follows, we will present results for different 
values of M and identify the optimum choice. 

As discussed above, the most relevant parameter of 
any trap model is the common barrier level b. For low 
energies this gives rise to the choice V(e) — b — e as 
the barrier height to leave a trap (i.e. MB) of energy e. 
To cope with the high-energy limit of V e //(e), discussed 
above, we choose V(e) — 1.0 for high e. The precise value 
is irrelevant for the results in the relevant temperature 
range. In total, we have V(e) = max(6 — e, 1). This 
function is sketched in Fig^Jx The average waiting time 
in a trap of energy e is given by 

r(e,T) = r exp(/3y( e )). (7) 

The value of To is a trivial scaling factor for the overall 
time scale and can be simply adjusted. 

C. Dynamic interaction between subsystems 

The ETM, introduced so far, starts from a collection 
of M totally independent subsystems. However, one 
should rather expect that a local relaxation process may 
somewhat influence the state of the adjacent subsystems. 
Thus, some kind of dynamic interaction should be in- 
cluded, i.e. a modification of the state of an adjacent 
subsystem as a result of a relaxation process. In real 
space this modification might correspond to a minor shift 
of particles in the adjacent subsystems. This shift may 
somewhat change the probability for a relaxation process 
in these subsystems. Actually, already in the original 
paper of Monthus and Bouchaud their master equation 
for the energy probability distribution has been supple- 
mented by an energy diffusion term which exactly takes 
into account the effect of relaxation processes on adjacent 
particles |30| . 

Here we introduce a simple ingredient of the ETM 
which reflects the relevant physics of the dynamic inter- 
action. In particular we will take care that the thermody- 
namics of the total system is not modified. This condition 
is non-trivial because any variation of states may directly 
influence the probability distribution of states and thus 
change quantities like the average energy. 

For the total system one may define the distribution 
(p(e, T) via 

^(e,T)cxp(e,r)/(r(e,r)). (8) 

It denotes the probability that after or before a transition 
a MB has the energy e. Thus, generation of MB accord- 
ing to the <p(e)-distribution (for convenience we suppress 



the dependence on T) yields the correct statistics of MB 
and thus the correct thermodynamics. Conceptually, the 
distribution <p(e) is close to the distribution of MB, i.e. 
G(e). Some minor (temperature dependent) deviations 
are present because in the present version of the trap 
model the relation V (e) = b — e is not valid for the high- 
energy traps with e > b — 1 . 

In what follows we generalize Eq|S]to the case M > 1. 
First we consider the case M = 2 but generalizations are 
straightforward. Let the subsystem 1 be in state e\ and 
the subsystem 2 in state e 2 . In the ETM the transition 
of the total system is achieved by a transition of one 
subensemble. Then the rate l/r(ei,e2) for a transition 
of the system is given by the sum of the individual rates, 
i.e. 

l/r(ei,e 2 ) = l/r( ei ) + l/r(e 2 ). (9) 

Furthermore due to the independence of the two sub- 
systems the combined Boltzmann probability p(e.\, e 2 ) is 
identical to the individual Boltzmann probabilities, i.e. 

p(ei,e 2 ) =p{ei)p(e 2 ). (10) 

Finally, the probability that after (or before) a tran- 
sition of the total system one finds energies e\ and e 2 is 
given by 

(/j(ei,e 2 ) oc p(ei,e 2 )/r(ei,e 2 ) 

oc p(ei)p(e 2 )(l/T(ei) + l/r(e 2 )) 

oc (/3(ei)r(ei)(/3(e 2 )r(e 2 )(l/r(ei) + l/r(e 2 )) 

oc (p(e 1 )ip(e 2 ){T(e 1 ) + T(e 2 )). (11) 

Due to the enormous dynamic heterogeneities (see be- 
low) the waiting times are distributed over many orders 
of magnitude. Thus, at a randomly chosen time one will 
typically find both subsystems in traps of very different 
waiting times T\ and r 2 . Let us assume that a jump pro- 
cess occurred in subsystem 1. Then it is very likely that 
the state before the jump was characterized by t± <C t 2 . 
Therefore one can write 

tp(ei,e 2 ) oc </j(ei V(e 2 )r 2 oc tp(ei)p(e 2 ). (12) 

Now we may define the dynamic interaction. It is mod- 
elled such that after a jump in one subsystem with prob- 
ability q a new trap is selected in the other subsystem. 
According to Eq^| the energy of the new trap in every 
subsystem has to be selected from the probability distri- 
bution p(e 2 ). This implies that the thermodynamics is 
not modified by the dynamic interaction. Furthermore, 
the principle of causality is directly implemented. Of 
course, in practice one would expect that every jump in- 
duces a minor rearrangement so that after 1/q jumps of 
the fast subsystem a total rearrangement with respect to 
the p-distribution has occurred. This gradual process, 
however, is very well approximated by the random pro- 
cess, introduced above. 
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Actually, for the low-energy limit, i.e. for r(e 2 ) oc 
exp(/3(6 — 62)), there exists a more intuitive rationaliza- 
tion of the factor p(e 2 ) in Eq^l It is reasonable to as- 
sume that the probability for subsystem 2 to change its 
energy from e 2 to e 2 will be proportional to t-pifi'i) and the 
Boltzmann factor exp(/?(e 2 — &•!))■ Therefore the proba- 
bility to be in state e 2 is proportional to </?(e 2 ) exp(/3(e 2 — 
e 2 ) tx <p(e 2 ) exp(-/5e 2 ) oc V5(e 2 )r(e 2 ) ocp(e 2 ). Of course, 
the original derivation of Eq^| is more general. 

For general M one obtains the relation 

ip(e lt ...,e M ) oc ip{e 1 )p(e 2 ) ■ ■ ■ p{e M ) (13) 

if a jump occurred in the first subsystem. Here we modify 
all subsystems m — 2,...,M with probability q accord- 
ing to the respective Boltzmann distribution p(e m ). Ev- 
idently, in reality the dynamic interaction between sub- 
systems may be somewhat more complicated. It will turn 
out, however, that this simple choice already captures the 
main effect of the interaction. 



D. How to compare the ETM dynamics with the 
BMLJ dynamics? 

So far we have introduced the ETM with three param- 
eters M, b, and q. Every state is characterized by the 
internal parameters e\, em- The time evolution of the 
system according to the ETM is generated in a straight- 
forward way. After a relaxation process in one of the M 
subsystems a new energy is selected for this subsystem 
according to the <p- distribution, defined in Eq|S] Further- 
more with probability q a new energy is selected for the 
other M — 1 subsystems according to the p — distribution, 
defined in Eq0J The total energy of all subsystems is 
denoted e l (the upper index counts the number of tran- 
sitions of the total system) . In the next step for every 
subsystem a waiting time T ac t(e m ) is calculated. It is 
drawn from an exponential distribution with the aver- 
age value r(e m ,T); see Eq[7| Then the subsystem with 
the shortest waiting time r act:m i„ is selected to perform 
a relaxation process and the time proceeds by exactly 
t % = T act , m i n . This procedure is then repeated and in 
the next step e l+1 and r l+1 are obtained etc. The out- 
put is a sequence of waiting times {t 1 } and corresponding 
energies {e 1 }. Exactly these two sequences can be also 
obtained from analysis of the equilibrium BMLJ dynam- 
ics. 

Both series together characterize in a detailed way how 
the system is exploring the PEL. They are to a large 
extent characterized by (i) the first and second moments 
of their respective distribution functions and (ii) by their 
correlation functions 

f,(n)^{^ a -^) 2 )/{^-^)) 2 ). (14) 

The variable n counts the number of MB transitions and 
/j 6 {t, e}. The denominator guarantees that / M (0) = 1. 
Basically, the f^,{n) express to which degree the waiting 



times and the MB energies are correlated after n MB 
transitions. Since the trap distribution in the ETM is 
chosen in order to reproduce the thermodynamics of the 
BMLJ system, the distribution of energies is naturally re- 
covered. Thus the observables (r) , (r 2 } , f T (n) and f e (n) 
may be taken for comparison. If the ETM is an appropri- 
ate model to characterize the PEL it should be possible 
to choose the model parameters M, b and q such that the 
BMLJ results of these observables are recovered for the 
relevant temperature regime between 0.45 and 0.6. 

It is important to clarify the influence of the dynamic 
interaction parameter q on the series of waiting times 
and energies, respectively. Since our procedure does not 
change the distribution (^(ei,e 2 ) (here for M=2) and 
since the waiting times are directly related to the en- 
ergies, the distributions of waiting times and energies do 
not depend on q. In contrast, due to the additional varia- 
tion of energies (and thus waiting times) one may expect 
that the decay of / M (n) becomes faster with increasing 
q. Thus, except for comparison of f^in) we may choose 
q = 0. 



IV. RESULTS OF THE MAPPING 

So far we have developed a minimum model for the 
PEL of a glass-forming system based on the trap model 
as the elementary unit. Now we compare the resulting 
ETM with the BMLJ system. As already mentioned this 
comparison involves the time series of waiting times as 
well as the time series of MB energies. 

A. Time series of waiting times 

First, we analyse the autocorrelation function / T (^)- 
The result for the BMLJ system is shown in FigOO Inter- 
estingly, it turns out that already after one step the new 
waiting time is basically uncorrelated to the old waiting 
time. This observation can be easily rationalized in terms 
of the ETM. From Eq. d (or Eq. U3 for general M) it 
follows that the new state of subsystem 1 will be drawn 
according to the distribution </?(ei). Since <p(e\) is peaked 
at much higher energies than p(e 2 ) it is very likely that 
also the next jump occurs in subsystem 1 (see below for a 
closer analysis of this effect). Since, however, two succes- 
sive energies in subsystem 1 are uncorrelated (by defini- 
tion) the same holds for the two successive waiting times 
for this subsystem. Since these waiting times determine 
the waiting time of the total system (because subsystem 
2 is basically immobile) one may conclude that in general 
two successive waiting times of the total system are to a 
very good approximation uncorrelated. Of course, it can 
be checked explicitly (data not shown). This non-trivial 
prediction fully agrees with the MD data for the BMLJ 
system. 

Now we discuss the first and the second moment of the 
waiting time distribution /(t) which, of course, strongly 
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Figure 5: The correlation function of waiting times in depen- 
dence of the number of MB transitions for the BMLJ system 
at T = 0.5. 

depends on temperature. As discussed above we may 
set q = for the comparison with the BMLJ system. 
In a first step we determine the respective values of b 
(for different values of M) from the condition that the 
temperature dependence of (t) should agree as well as 
possible with the BMLJ data. 

There are two different ways to estimate the value of 
(t). On the one hand, one can calculate the average of all 
waiting times, observed during a simulation run of time 
tsim ■ On the other hand, one may determine the number 
of transitions n and take t B i m /n as the average waiting 
time. As shown in Appendix I the second choice can be 
applied for any simulation time i s j m . In contrast, the first 
choice is only applicable for very large t S i m . Therefore 
we prefer to use the second choice. 

In FigElwe show the temperature dependence of the 
average waiting time (r(T)) as obtained from the MD 
simulations. These data are used to determine an opti- 
mum value of the common barrier level b for the different 
choices of M. The resulting predictions for the ETM are 
also included in Fig|5] The values for 6, obtained from 
this fitting, are listed in Tab.l. The agreement with the 
BMLJ data is not perfect but this single fitting parame- 
ter b is enough to reproduce the temperature dependence 
in the relevant temperature range between 0.45 and 0.6 
in a semi-quantitative way. By increasing the value of 
b it would have been possible to improve the agreement 
for low temperatures whereas the deviations at higher 
temperatures would have become stronger. In any event, 
the temperature dependence of the ETM for the differ- 
ent values of M looks very similar. Deviations at higher 
temperatures are expected because in this regime flow- 
like processes become more important which, of course, 
cannot be fully captured by the ETM. 

So far we have analyzed the average waiting time. At 
a given temperature it has contributions from MB with 
very different energies. Therefore it is more informative 
to analyze the waiting time in dependence of temperature 
and energy. As expressed by Eq. [S] there is a significant 
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Figure 6: The average waiting time of MD simulations as 
compared with the predictions of the ETM with different val- 
ues of M. The resulting values of 6 and to for the different 
M are given in Tab.l. 



M 


b 


To 


2 


-143.1 


160 


3 


-93.3 


43 


4 


-68.3 


14 


6 


-43.4 


4 



Table I: The values of the barrier level b and the prefactor to, 
obtained from the fitting as discussed in Fig[S] 



dependence of the average waiting time on the MB en- 
ergy in the BMLJ system. For this purpose we consider 
the energy-dependence of the quantity Tln(r(e, T)). For 
a simple trap model (see Eq. [7J it is linear with slope -1. 
What is the result for the ETM? The result for M = 3 is 
shown in FigEK for the energy interval [-300,-295]. Inter- 
estingly, for all temperatures (T=0.45;T=0.5;T=0.6) a 
constant slope is observed which furthermore is identical 
for the three curves. For the energy interval [—300, —295] 
the data for the different temperatures can thus be con- 
sistently described by 

T ln(r(e, T)) « const{T) - ce (15) 

with c = 0.33. We have repeated the same analysis for 
the different values of M. One observes the same behav- 
ior except for different values of c; see Tab. 2. 

Actually, the value of c can be also estimated from 
general arguments which have been already partly pre- 
sented in pol ] . For M subsystems and individual energies 
(ex, cm) in the limit of low temperatures the activation 
energy is given by min[T/(ei), ...V(e M )] = V(e max ) where 
e-max is the largest value of the {e m }. To calculate the 
apparent activation energy V e f / (e) for fixed total energy 
e, one has to average this expression over different real- 
izations of the e m under the constraint e = e% + ... + eu- 
As a first approximation, however, one may take e max « 
e/M + Ae. The second term accounts for the deviation 
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M 


c 


2 


0.44 


3 


0.33 


4 


0.25 


6 


0.17 



Table II: The slope c obtained from analyzing the energy 
dependence of Tln(r(e,T)> in Figd 



°f e max from the average value e/M. It can be expected 
that Ae only weakly depends on e. Then one has 

Veff{e) = V(e max ) = V(e/M + Ae) « const - (1/M)e 

(16) 

Thus one expects that the slope c is close to 1 /M which 
is in excellent agreement with the numerical results in 
Tab.2. 

In analogy to the ETM we have calculated the quantity 
Tln(r(e, T)) for the BMLJ system for the same tempera- 
tures. The results are displayed in Fig03. Interestingly, 
the data also display a linear energy-dependence. There- 
fore it is possible to directly extract the value of M if the 
BMLJ data are interpreted in terms of a ETM. For the 
two lowest temperatures the resulting slope is closest to 
M = 3. For reasons of simplicity we also choose M = 3 
for the T = 0.6. The weak temperature-dependence of 
the slope will be further discussed below. We note in 
passing that for a monatomic Lennard-Jones system of 
32 particles the average barrier height linearly depends 
on energy with slope one. This would correspond to M=l 
and is thus qualitatively compatible with the present re- 
sult for a larger (and somewhat different) system |46l |. 

In the next step we want to elucidate the nature of the 
fluctuations. A natural choice would be the quantity 

S ( T ) s - (r)2 ^ (17) 

including the first and the second moment of the waiting 
time distribution f(r). It can be regarded as a measure 
for the normalized variance of the waiting time distri- 
bution. Around T c one approximately finds /(r) oc t~ 2 
for large r. Without a cutoff of /(r) at long times the 
variance would diverge for this distribution. This implies 
that the variance is extremely sensitive on the cutoff- 
behavior of /(t). For temperatures around T c this cutoff 
can be predicted from Eq[7|by setting e = —306 as an 
estimate for the low-energy cutoff of the energy land- 
scape. It turns out T rnax w 10 5 (r) which corresponds to 
roughly 10 10 MD steps. Only for simulation times much 
longer than T max and thus many orders of magnitude 
longer than possible by present computer technology the 
variance of the observed waiting time distribution can 
be determined from simulations around T c . This will be 
explicitly demonstrated in the next Section. 

Evidently, the very long waiting times result from MB 
with very low energies. Therefore it may be very instruc- 
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Figure 7: The quantity Tln(r(e,T)) for different energies 
and temperatures in the energy range [-300,-295]. a) Results 
for the ETM for the example M = 3. The resulting slopes 
for all relevant values of M are listed in Tab.2. b) Results for 
the BMLJ. The resulting temperature-dependent slopes are 
included in the Figure. 



tive to analyze the fluctuations of waiting times of MB, 
restricted to energy e. This restricted quantity is de- 
noted as S(e,T). For very low energies one may expect 
that due to the extremely small number of contribut- 
ing MB the estimated variance will be smaller than the 
true variance. Only for higher energies S(e,T) can be 
estimated without systematic errors from finite simula- 
tion times. We restrict ourselves again to the interval 
[—300, —295]. S(e,T) is expected to show a sensitive de- 
pendence on the nature of the ETM. If, for example, the 
total BMLJ system could be described as a single trap 
system (M = 1), MB with energy e would be described 
by a single relaxation rate, yielding 5(e) = 1. For M > 1 
the total energy can be decomposed in different ways, re- 
sulting in different values of r(ei, ...,ej\/) for fixed total 
energy e% + ... + ejvr; see Eq|5J The resulting distribution 
is a sum of different exponentially distributed functions, 
giving rise to 5(e) > 1. The data for the BMLJ sys- 
tem for T = 0.45; 0.5; 0.6 are shown in the inset of FigEl 
for the interval e G [—300,-295]. The average value of 
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Figure 8: Inset: Energy-dependence of the normalized vari- 
ance S(e, T) for different temperatures for the BMLJ sys- 
tem. Main panel: Comparison of the respective average values 
S {T) with the predictions of the ETM for M = 3 and M = 6. 



S(e,T) in this interval is denoted Sq(T). It is plotted in 
Fig|S]for the three temperatures. Furthermore we have 
added the predictions for S {T) of the ETM for M = 3 
and M = 6. It turns out that the data are inconsistent 
with the predictions of the ETM for M — 3. Rather a 
good agreement with M = 6 is suggested. As will be 
discussed further below a very natural mechanism may 
exist which reduces the size of the fluctuations of waiting 
times for a real system as compared to the predictions of 
the ETM. Therefore we propose that M ss 3 as obtained 
from the analysis of (r(e, T)) is indeed the relevant value 
for the ETM. 



B. Sequence of energies 

In a last step we analyse the sequence of MB ener- 
gies. Since the energy distribution of the real system is 
naturally recovered by the trap model (see above) the 
non-trivial aspect is related to possible correlations of 
successive MB energies as measured by the autocorre- 
lation function f e (ji). The function f e (n), obtained for 
the BMLJ system, is shown in Fig|5J Since for T = 0.5 
one has T a « 20 (t) (r a determined from the incoherent 
scattering function) the energy correlation decays on the 
time scale of the a-relaxation. 

Dramatic differences are visible when comparing f e (n) 
with the predictions of the ETM (M=3)without including 
the dynamic interaction, i.e. q — 0. Whereas the energies 
for the BMLJ system significantly decorrelate after 40 
steps, basically no decorrelation is seen for the ETM. 
This can be easily rationalized. In the typical scenario, 
discussed above, very often there exists one subsystem 
which is much faster than the other two subsystems and, 
correspondingly, has a relatively high energy. This fast 
subsystem jumps and, on average, will acquire a new 
waiting time which is still much faster than the waiting 
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Figure 9: Decay of the energy-autocorrelation function in de- 
pendence of the number of MB transitions at T=0.5. Shown 
are the MD data for the BMLJ system (circles) and the predic- 
tions with (q = 0.12) and without (q = 0.0) dynamic coupling 
for the case M = 3. 



time of the other two subsystems. Thus the energy of the 
slower systems will remain unchanged even after a large 
number of transitions. This automatically implies strong 
correlations for successive total energies of the ETM. 

In contrast, a finite dynamic interaction results in 
additional reorganisation processes of the slow subsys- 
tems, thereby changing their energy according to the p- 
distribution. Although p(e, T) is shifted to lower ener- 
gies as compared to <^(e, T) there will be the chance to 
acquire a relatively high energy after a few transitions of 
the fast subsystem and start to contribute to the relax- 
ation. Therefore one can expect that a finite q will give 
rise to a considerable decorrelation of energy. This can 
be seen in FigEl where data for q = 0.12 are included. 
With this value the initial part of the decorrelation can 
be reproduced. Taking the strong dependence of f e (n) on 
the value of q into account one may estimate the uncer- 
tainty of this value by ±0.02. Actually, a similar analysis 
for the two other temperatures yield q- values in the same 
interval although the data suggest a minor temperature- 
dependence (q(T = 0.6) « 0.14; q(T = 0.45) « 0.11). 

Of course, this mechanism of dynamic interaction is far 
too simple to explain the complex coupling between dif- 
ferent relaxation modes in supercooled liquids. In partic- 
ular, it turns out that finally the decorrelation, generated 
by the dynamic interaction, is too strong. From a physi- 
cal point of view one would expect that states with lower 
energies are less sensitive to fluctuations in the neighbor- 
hood than states with somewhat higher energies. This 
implies a dependence of the strength of the dynamic in- 
teraction on energy. Qualitatively, this property would 
explain the remaining deviations in FigEI Unfortunately, 
when introducing an energy-dependent q- value one would 
be forced to modify the simple rule based on Ea ll3l in or- 
der to obtain the same thermodynamic properties. This 
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is beyond the scope of the present work and would not 
yield any relevant new insight into the nature of the PEL 
of the BMLJ system. 



V. DISCUSSION 

A. Length scales and finite-size effects 

Analysing the time series of waiting times the BMLJ 
system with 65 particles behaves to a first approxima- 
tion as if it were a superposition of M « 3 independent 
subsystems. Each subsystem contains on average 20-25 
particles. However, performing MD simulations with this 
very small number of particles (and periodic boundary 
conditions) one would observe strong finite size effects. 
Thus, it is not possible to perform simulations on the 
level of the elementary units of the ETM. In this sense 
it is impossible to see the elementary trap system di- 
rectly. Rather the properties of the individual subsys- 
tems, mainly expressed by the value of the barrier 6, 
have to be extracted from the superposition of at least 
2-3 subsystems, corresponding to roughly 65 particles in 
real space. Why does a PEL of a BMLJ system with only 
20 particles show significant finite size effects? First, the 
presence of dynamic interaction directly shows that the 
coupling among adjacent regions of the supercooled liq- 
uid is essential for the dynamics. Second, one can imagine 
that the real-space pattern of those particles which move 
together during one MB transition may be shaped in an 
string- like manner [i^, [5(| • Therefore the minimum size 
system without major finite size effects has to be larger 
than the elementary subsystem of the ETM. Otherwise 
the real-space patterns cannot be formed appropriately. 

Of course, one could repeat the analysis also for larger 
systems. Evidently, larger values of M would have been 
found. In any event, since a BMLJ system with N = 130 
particles behaves like a superposition of two systems with 
N = 65 particles one would end up with a ETM which 
is a trivial extension of the present system. 



B. How to rationalize the size of waiting time 
fluctuations at constant energy? 

The only observable which did not agree with the pre- 
dictions of the ETM was the second moment of the wait- 
ing time distribution. The fluctuations at constant en- 
ergy as expressed by So (T) are smaller than expected for 
M = 3. Is there a simple physical picture which may 
reconcile these results? 

Recently a detailed real-space analysis of the nature 
of the MB transition has been performed for the BMLJ 
system |39j . In particular the number of particles, par- 
ticipating at such a transition, has been studied. This 




Figure 10: Sketch of the size fluctuations of a ETM with 
M = 3. 

number can be quantified by the participation number 

Here the sum is over all large particles of the binary 
Lennard-Jones system, dri is the displacement of the 
i-th particle and dR 2 is the total squared displacement 
of all large particles. In case that only one particle is 
moving one has z — 1, in case of identical movement of 
n particles one gets z = n. It turns out that for different 
MB transitions the participation number z can have very 
different values. The resulting distribution p(z) has sig- 
nificant contributions between z w 8 and z ~ 25. Since 
the absolute numbers somewhat depend on the definition 
of z |3^| we mainly stress the width of p(z) rather than 
the absolute numbers. 

In the ETM, as introduced above, we have introduced 
M identical subsystems. If this corresponded to the real 
world of BMLJ systems, one would expect a very narrow 
distribution p(z). The width of p(z) therefore implies 
that the assumption of identical subsystems cannot hold. 
In any event, it is natural to assume that at a given time 
the set of subsystems displays different sizes; see FigUHl 
for a simple sketch. 

To understand the effect of these fluctuations on the 
waiting time distribution one may consider the case of 
M identical subsystems for fixed total energy e. The 
maximum waiting time of the total system is realized if 
by chance all M subsystems have the energy e/M and 
the corresponding waiting time r e . Otherwise there is 
at least one subsystem with higher energy and (on aver- 
age) shorter waiting time as compared to r e . Thus the 
presence of the long-time wing of the waiting time dis- 
tribution at constant energy is critically related to the 
presence of subsystems with identical waiting times. Al- 
lowing for size fluctuations it is likely that one subsystem 
is smaller than the other. Since on average smaller sys- 
tems are faster than larger systems (see, e.g., Tab. 2) one 
can expect that for a given total energy e this subsys- 
tem and thus the total system has a waiting time much 
smaller than r e . Therefore size fluctuations strongly re- 
duce the probability to have a situation with a very long 
waiting time for the whole system. As a consequence, 
these fluctuations strongly suppress the long-time tail of 
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the waiting time distribution, thereby reducing higher 
moments of this distribution. One may speculate that 
this is the reason that the fluctuations for the BMLJ are 
somewhat smaller than expected from the average wait- 
ing times. Of course, additionally also the value of M 
may somewhat differ for different MB. 

C. Implications for dynamic heterogeneities 

The presence of short and long waiting times can be di- 
rectly related to the concept of dynamic heterogeneities, 
observed in many experiments. As discussed above, a 
relevant measure for the presence of dynamic hetero- 
geneities is the value of S(T). On the level of the elemen- 
tary subsystem in the ETM the dynamic heterogeneity is 
exclusively related to the presence of different trap depths 
which, via Eq[7| gives rise to a broad distribution of wait- 
ing times. For the BMLJ system with 65 or more parti- 
cles, i.e. (in the language of the ETM) for a superposi- 
tion of a few elementary subsystems, the dynamic hetero- 
geneities can be formally divided into two parts. First, al- 
ready at constant total energy e one has a range of differ- 
ent waiting times as reflected by Sq(T) > 1 (for M > 1). 
Second, the average waiting time (r(e, T)) depends on en- 
ergy. Comparing the size of ^/ Sq(T = 0.45) « 3 with the 
range of the energy dependence of (r(e, T)) it becomes 
obvious that for 65 particles the dynamic heterogeneities 
are mainly determined by the energy dependence of the 
average waiting time. 

It is possible to characterize the properties of dynamic 
heterogeneities somewhat closer. As discussed before, the 
dynamic heterogeneities or, equivalently, the width of the 
waiting time distribution is captured by the normalized 
variance S(T) (see Eg 1170 . In the limit of very long sim- 
ulation times t s im as compared to the range of waiting 
times, the equilibrium value of S(T) can be extracted 
from the simulations. How to find an appropriate way 
of extracting corresponding information at shorter simu- 
lation times? As already discussed in the context of the 
determination of (r(T)) it is convenient to analyse the 
number of MB transitions. Beyond the determination of 
the average number of MB transitions one can estimate 
its variance ((n — (n)) 2 ) by comparing the number of MB 
transitions during a fixed time t S i m of several indepen- 
dent simulations. Then one may define 

S{TM m )= {{n - { $ ?) t -^. (19) 

As shown in Appendix II in the limit of long t s \ m the 
quantity S(T,t s i m ) approaches the equilibrium value of 
S(T). For shorter times S(T,t S im) reflects the normal- 
ized variance of that part of the waiting time distribution 
which can be accessed on the time scale of the simula- 
tion tsim- Thus any dependence of S{T,t S i m ) on t a i m 
indicates that the finite simulation time introduces an 
artificial cutoff for the waiting time distribution. As 
shown in FigEl for T < 0.5 the value of S(T,t sim ) 
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Figure 11: The dependence of the normalized variance 
S(T,t s im) on simulation time t s i m and temperature for the 
BMLJ system. The straight lines are the estimates of the 
ETM for the respective limiting values. The estimates for the 
ETM have been scaled by a factor, reflecting the deviations 
for Sb(T) between the ETM and the BMLJ system. 



for the BMLJ system still shows a strong increase for 
tsim ~ 1000 (t) ~ 20t q which already corresponds to rel- 
atively long simulation times. 

It may instructive to compare the long-time limit of 
S(T,t s im) for the BMLJ system with the predictions of 
the ETM. We have estimated the latter value by taking 
(using M = 3) S(T) ■ [S ,bmlj(T)/ S ,etm{T)}. The 
second factor corrects for the fact that the actual dy- 
namic heterogeneities of traps with identical energy are 
smaller than the predictions of the ETM; see FigEl The 
results are included in Fig^] For the highest temper- 
ature T = 0.6 the long-time limit of S(T,t S i m ) can be 
reached for t S i m w 1000(r). Good agreement with the 
prediction of the ETM is observed. The large estimates 
of S(T) within the ETM for the two lower temperatures 
reflect the broad dynamic heterogeneity at low tempera- 
tures. They only fully show up in the limit of very long 
simulation times. 

This result has interesting consequences for the plan- 
ning of simulations. The value of S(T, t s i m ) is a measure 
how precisely the true value of (t) can be extracted from 
the simulation. More specifically, from L independent 
runs the value of (r) can be extracted with a precision 
of y/ (S(T, t sim ))/(Lt sim / (r)). Here we have assumed 
that a simulation of length t S i m implies the drawing of 
tsim/{ T ) independent waiting times. If S(T,t S i m ) does 
not depend on t S i m one could either perform a long sim- 
ulation or several short simulations to obtain the same 
quality of (r). Due to the time-dependence of S(T,t S i m ) 
it is more favorable to choose t S i m relatively short and 
thus perform many independent runs for fixed total com- 
puter time. Alternatively, one can also perform a few 
simulations for a correspondingly larger system and ob- 
tain the same quality for the average waiting time. Since 
the diffusion constant is directly related to (r) the same 
conclusions hold for the determination of the diffusion 
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constant. 

Recently, the variance of the mean square displacement 



((r 2 



)))/(') nas been obtained from comparison 



of independent simulation runs 23]. On a superficial 
level this ratio resembles the quantity ((n — (n)) 2 )/(n) 2 
analyzed above. A possible mapping between the ob- 
servables r 2 and n is further suggested by the validity of 
Eq. ^ Closer analysis, however, reveals that this map- 
ping is not possible. For a simple random- walker in ID 
one has p(r 2 )dr 2 = exp(— r 2 / (r 2 ))dr 2 [2^. In contrast, 
according to Appendix II the the value of n displays a 
gaussian distribution around (n) . As discussed in [23 the 
information content of the variance of the mean square 
displacement is more related to the size of dynamic het- 
erogeneities whereas in the present case we are sensitive 
to the distribution of waiting times, i.e. to the degree of 
dynamic heterogeneities. In any event, the common idea 
is to use independent simulations to extract important 
information about the non-trivial dynamic properties of 
supercooled liquids. 

Following the general ideas, underlying the Adam- 
Gibbs scenario or alternative descriptions of the physics 
of supercooled liquids at low temperatures |§3> 13 EH 
one might expect a growing length scale of the cooper- 
ativity range. In the present approach the length scale 
is reflected by the slope of Tln(r(e,T)). Indeed, some 
temperature dependence has been observed, changing the 
slope by roughly 30 % when changing the temperature 
from T c (0.45) to 1.33T C (0.6). In terms of the ETM this 
implies a weak temperature dependence of the coopera- 
tivity size. Actually, the slopes for T = 0.45 and T = 0.5 
were compatible with M = 3. 

It may be instructive to compare this weak temper- 
ature dependence with that of other observables which 
characterize length scales of dynamic heterogeneities. 
Analogous behavior is observed for the average value of 
the participation number z; see Eg 1181 This value is ba- 
sically constant when comparing T — 0.5 and T = 0.435 
|39|. In contrast, a strong temperature-dependence is 
observed for quantities which are determined from the 
space-time behavior of particles. For example, the maxi- 
mum of the generalized susceptibility , characterizing the 
size of cooperativity regions, varies by roughly a factor 
of 5 in the temperature regime analyzed in this work; see 
e.g. the behavior of \m hi Rcf. 56] . How to rationalize 
this apparent discrepancy? In the language of MB, dy- 
namic heterogeneities have different facets: (i) The spa- 
tial range of particles moving during individual MB tran- 
sitions, (ii) Spatial correlations of subsequent transitions 
|39j . Both aspects are part of the ETM. The first aspect 
is directly reflected by the value of M. The second aspect 
is related to the question how many successive hops are 
on average performed by the same subsystem until a dif- 
ferent subsystem relaxes. This value is denoted n succ . It 
can be easily calculated from simulation of the ETM. As 
shown in Fie ll 21 n is strongly temperature-dependent 
and increases within the temperature regime of interest 
by a factor of 6. In the space-time analysis one is sensi- 
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Figure 12: The average number n succ of successive hopping 
processes of the same subensemble in the ETM with M = 3 
for different temperatures. 



tive to the sequence of fast processes, i.e. the sequence of 
MB transitions. Thus one may expect that the increase 
of the generalized susceptibility with decreasing temper- 
ature is at least partly due to the increase of n succ . On 
a qualitative level similar effects have been already re- 
ported in |39| where the formation of macro-strings out 
of micro-strings has been observed. 



VI. SUMMARY 

Many observations for the BMLJ system are consis- 
tent with the notion that the PEL can be characterized 
as a superposition of individual trap-like systems where 
the individual traps are identified as MB (and not as in- 
herent structures). The main observations are (i) the 
random-walk nature of the dynamics in phase space, (ii) 
the exponential waiting time distribution to leave MB, 
(iii) the temperature independent distance of subsequent 
MB, (iv) the absence of correlations of subsequent MB 
waiting times, and (v) the linear energy-dependence of 
Tln(r(e, T)) for all temperatures. Analysis of the slightly 
temperature-dependent slope of Tln(r(e,T)) yields in- 
formation about the number M of subsystems and the 
common barrier height b. 

From the predictions of this simple model two devia- 
tions from the behavior of the BMLJ system have been 
observed. First, the energy correlation function decays 
faster than expected from the superposition of traps. 
This observation directly indicates the presence of some 
interaction between these subsystems, expressed by the 
dynamic interaction parameter q, thereby completing our 
definition of the extended trap model (ETM). Second, 
the dynamic heterogeneities of MB with constant e were 
smaller than expected. This may be a natural conse- 
quence of size fluctuations of the subsystems. In principle 
it would be possible to include size fluctuations into the 
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ETM by introducing a new parameter. This is, however, 
beyond the scope of the present work. 

Our results may be helpful in several ways. The com- 
parison of the BMLJ dynamics with the predictions of 
a superposition of appropriately chosen trap models may 
elucidate important properties of the PEL of a prototype 
glass-forming system. Furthermore, this analysis may 
serve as a fresh look onto the properties of dynamic het- 
erogeneities. Finally, the ETM has the ability to repro- 
duce also non-equilibrium properties at least on a semi- 
quantitative basis. For example it may be feasible to 
reproduce complex aging experiments and to obtain a 
simple physical picture of the resulting observations. For 
this purpose it is helpful that the relevant length scale for 
the ETM, as expressed by the value of M, is only weakly 
temperature dependent. 

Appendix I 

We consider a simple system which can be in differ- 
ent states with weights tpi. The relaxation process in 
each state is characterized by an exponential waiting time 
distribution with average r^. Thus the resulting average 
waiting time is given by (r) = 'Y^i L Pi T i- I n equilibrium 
the probability to be in state i is given by pi — (fiTi/ '(r). 
First, we consider a simulation of time t S i m which is 
much smaller than any of the Tj. In particular one has 
tsim <C (t) . The goal is to estimate the value of (r) from 
the outcome of the simulation. As described in the main 
text we may proceed in two different ways. One option 
is to determine the average waiting time. Since t S i m is 
very small at most a single relaxation process will be ob- 
served. For calculating the average waiting time (in the 
ensemble average) the most convenient way would be to 
restrict oneself to those instances where at least a single 
relaxation process has occurred. Due to the smallness of 
t S i m it may occur at any time with equal probability so 
that the average waiting time will be t S i m /2. This value 
is, of course, much smaller than the true average waiting 
time (r). Actually, only if t S i m is much larger than all 
Tj a consistent estimate is possible. The second option is 
to determine the average number (n) of relaxation pro- 
cesses. Weighting the initial state with its probability 
Pi, i.e. assuming equilibrium conditions, and taking into 
account that the chance of a relaxation process during 
time t S im is given by t S i m /Ti one ends up with 

/ v X > t s im \ ^ tsim t s im 

(n)=2^Pi—=2^^i 1 -r= 1 -r- (20) 



Thus t s im I '(n) is an exact estimate of the average waiting 
time (t) . This result also holds for larger values of t S i m 
because the large time interval can be decomposed into 
small intervals and the average value of (n) will simply 
scale with the number of small time intervals such that 
(n) = t s im I '{t) will hold for arbitrary f s j m . 
Appendix II 

Of conceptual interest is the variance a 2 of the wait- 
ing time distribution. We show that it is related to the 
variance er^, obtained from a set of independent runs. 
According to the central limit theorem the probability 
that it takes a time t to perform n jumps is given by 

pft|»)ocex P ( ~ (t 2 ~y )2 ). (21) 

Here we used that successive jumps are independent from 
each other so that we can consider a drawing of n in- 
dependent elements from the waiting time distribution. 
Fixing the value of t as t S i m this probability can be rein- 
terpreted as the probability p(n\t) that exactly n jumps 
occur during time t S i m . Thus one obtains 

p{n\t sm ) oc exp ^ 2na 2 /{T) 2 J ■ ( 22 ) 

In the limit of large t S i m the n-dependence of the normal- 
ization factor can be neglected. Furthermore the factor n 
in the denominator of the exponential can be substituted 
by (n) — tsim/ (t). Therefore the variance of n, i.e. <r n , 
can be written as 

(n) 2 (r)(n) t sim <r} 2 1 ' 

This directly leads to EaU^l 
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